AD-A 142  894  CONVERGENCE  OF  BIVARIATE  CARDINAL  INTERPOLATIONS)  1/1 

Wisconsin  univ-madison  mathematics  research  center 

C  D  BOOR  ET  AL.  MAV  84  MRC-TSR-2686  DAAG29-80-C-0041 

F/G  12/1 


UNCLASSIFIED 


NL 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANOAROS-1963-A 


UNIVERSITY  OF  WISCONSIN-MADISON 
MATHEMATICS  RESEARCH  CENTER 

CONVERGENCE  OF  BIVARIATE  CARDINAL  INTERPOLATION 

1  1  2 
Carl  de  Boor  ,  Klaus  Hollig  and  Sherman  Riemenschneider 

Technical  Summary  Report  #2686 
May  1984 

ABSTRACT 

He  give  necessary  and  sufficient  conditions  for  the  convergence  of 
cardinal  interpolation  with  bivariate  box  splines  as  the  degree  tends  to 
infinity. 


AMS  (M06)  Subject  Classifications:  41A05,  41A15,  41A63 

Key  Words:  bivariate,  box-splines,  interpolation,  convergence,  exponential 
type,  Whittaker  operator 

Work  Uhit  Number  3  (Numerical  Analysis  and  Scientific  Computing) 


^Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-80-C-0041. 
^ supported  by  NSERC  Canada  through  Grant  #A  7687. 


SIGNIFICANCE  AND  EXPLANATION 

This  is  a  follow-up  on  thfr  MRC  TSR  #2485  in  which  we  introduced  and 
studied  interpolation  by  a  linear  combination  of  translates  of  a  bivariate  box 
spline  on  a  three-direction  mesh.  This  is  of  interest  because  these  box 
splines  are  not  just  tensor  products  of  univariate  B-splines  but  are  genuinely 
bivariate/  yet  are  true  generalizations  of  the  univariate  cardinal  B-spline. 
This  allows  one  to  be  guided  by  Schoenberg's  highly  successful  analysis  of 
univariate  cardinal  splines/  while  at  the  same  time  struggling  with  a  more 
complicated  setup. 

The  specific  task  of  the  present  report  is  the  derivation  of  necessary 
and  of  sufficient  conditions  for  the  convergence  of  the  box  spline 
interpolants  as  the  degree  goes  to  infinity.  The  conditions  are  stated  in 
terms  of  the  Fourier  transform  of  the  interpoland. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 
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Denote  by  n*  the  "middle”  component  of  n,  i.e.,  the  second  number 


in  any  ordering  of  n^,  n2»  n^.  We  write 

n  ♦  N 


if  a  sequence  n(m),  ra  e  W,  satisfies 

(nl)  n'(ni)  ♦  *  as  m  ♦  «, 

(n2)  lim  -  N  e  [0,-J3. 
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Further,  we  assume  that 

(n3)  j n |  n^  +  n2  +  n3  <  c(n’)c 

where  c  is  some  positive  constant*  Examples  of  admissible  sequences  are 

n(m)  -  (m,2m,3m)  with  N  -  (1/2,  1,  3/2)  , 
n(m)  *  (1,  m,  m^)  with  N  «  (0,1,*). 


The  assumption  (n3)  excludes  degenerate  cases  like  n(m)  -  (l,m,m!). 

The  role  of  the  Interval  (-*,*)  is  played  by  certain  domains  fljj 

3 

corresponding  to  the  limit  of  the  sequence  n.  For  N  e  [0,«)  they  are 
defined  by 


(4)  iljj  {2xx:  0  <  aN  ^(x)  <  1  for  j  e  j) 

where  J »■  {±(1,0),  ±(0,1),  ±(1,-1)}  and  for  x  -  (u,v),  j  -  (k,i). 


(5) 
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u+v+k+Jt 


Clearly,  the  set  is  bounded  by  the  curves  TN  ^  :•  (2xx:  aN  j(x)  *  l}, 


j  e  J.  If  one  of  the  components  of  N  equals  <*>  the  sets  as  well 

as  the  curves  rM  have  to  be  Interpreted  as  the  appropriate  limits 
J 

(cf.  Proposition  2).  A  qualitatively  correct  picture  of  Is  given  in 

Figure  1.  Figure  2  shows  a  few  special  cases.  Of  particular  Interest  Is 
the  symmetric  case  N  -  (1,1,1). 
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Figure  2 


A  detailed  discussion  of  the  properties  of  the  sets  is  given  in  (3] 

We  merely  note  that  they  are  fundamental  domains,  i.e.,  up  to  a  set  of 
measure  zero,  their  translates  2vj  +  0^,  j  el2,  form  a  partition  of  R2 


Our  first  result  is  an  extension  of  Theorem  5.2  in  (3)  to  include 


interpolation  of  data  with  power  growth  as  was  done  in  [7)  for  the 
univariate  case. 

Theorem  1.  Assume  that  the  Fourier  transform  of  f  is  a  tempered 
distribution  with  supp  {  c  If  the  sequence  n  satisfies  (nl)  -  (n3), 

then,  for  any  a  the  partial  derivative  DaInf  of  the  cardinal 

2 

lnterpolant  converges  locally  uniformly  on  I T  to  Daf,  as  n  ♦  N. 

As  for  the  univariate  case,  the  converse  of  the  above  Theorem  holds 
with  "c  replaced  by  "c 

l 

Theorem  2.  Assume  that  the  segugfence  n  satisfies  (nl)  -  (n3).  If  a 
sequence  of  cardinal  splines  an  e  Sn  converges  locally  uniformly  to  f 
and  if  |*n(x)|  £  c(l  +  |xj ]C  for  all  n  and  some  c  >  0,  then 
supp  { 

He  may  relax  the  assumption  (n2).  Clearly  any  subsequence  of  n 

also  satisfies  (nl)  and  (n3).  If  {NaJ  are  the  limit  points  of  the 

sequence  n/n*  then  one  has  to  replace  the  set  in  the  Theorems  by 

n  a.  .  The  figure  below  shows  the  intersection  ft  and  the  union  SI  of 
«  a 

all  possible  limit  sets. 


I iRure  3 


2.  Proofs.  We  assume  throughouc  that  the  sequence  n  satisfies 
(nl)  -  (n3).  By  c  we  denote  various  positive  generic  constants  which 

do  not  depend  on  n.  These  constants  may  change  even  within  the  same 
line.  Further,  we  set 


dn(x)  dist  (x ,  5>Qn) 


and  denote  by  ^  the  characteristic  function  of  the  set 

Denote  by  e  Sn  the  fundamental  spline  which  Interpolates  the 
2 

data  6^  k  e  TL  .  It  is  easily  seen  [3]  that  I,n  decays  exponentially 
at  infinity.  Therefore,  if  we,  e.g.,  assume  that 

(2)  jf(x)j  <  c(l  +  |x|)c  , 

then  we  can  write  the  cardinal  lnterpolant  in  Lagrange  form 


If  «  l  ,  f(j)Ln(-J). 


The  proof  of  Theorems  1  and  2  is  based  on  the  following  estimate  for  the 

Fourier  transform  of  L„  which  will  be  derived  at  the  end  of  this 

n 

section. 


Theorem  3.  For  any  e  >  0  and  a  e  l*  there  exists  n^  such  that  for 
n'  £  n^  and  d(x)  >  e 


jDa(Ln(x)  -  X^x))]  <,  (l  +  cd^x))”11  . 


Proof  of  Theorem  1.  Denote  by  S  (4]  the  space  of  rapidly  decreasing  test 

2  a  a 

functions  +  elR  •  The  assumption  f  e  S’  and  supp  f  c  0^  implies  (2) 
and  hence  the  representation  (3)  is  valid  for  the  cardinal  lnterpolant. 
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Since  is  a  fundamental  domain  which  contains  supp  f ,  the  values 

f(-j)  -  (2iO"2<£,  e^*> 


are  the  Fourier  coefficients  of  £,  Therefore  £ R  converges  in  S'  to 
the  periodic  extension  of  £: 

f°  I,  f ( •  +  2»J). 

it*  , 

This  means  that  there  exists  y  eS+  such  that  for  any  ♦  e  S, 
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where 


<f  -  fK.  *>|  -  0(1)  •  1*1^,  as  K 


l*lv  max  sup,  ly^Ky)! 
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Note  that  (5)  implies 


(5')  |<fK,*>|  <  cl*lY, 

uniformly  in  K.  Putting  f(y)  :*  (iy)°eixy  we  can  mite  the  difference 

(2»)2  D°(f  -  Inf)(x) 

in  the  form 
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<£,  ( l-£  ) ♦>  -  lim  <#„-#,  £  ♦>. 
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Since  supp  £  c  the  first  term  can  be  estimated  using  Theorem  3 


with  e  dist(supp  i,  dfi^)/2.  For  the  second  term,  we  choose  a 
cut-off  function  w  e  S  with  supp  w  :»  U  {y:  d^(y)  <  e})  and 

with  w(y)  •  1  for  y  e  Since  ((2nj)  +  supp  £)  n  supp  w  ■  0  for 

j  *  0  and  supp  f  fl  supp(l-w)  *  0  we  have 

<i  -  /K,  Ln^>  -  <f“#  -  fK,  wLn^>  +  <-fR,  (1-w )tj>. 

The  first  term  tends  to  zero  as  K  -*•  «.  As  to  the  second  term  note  that 
dist  (supp(l-w),  fijj)  e/2,  which  by  (4)  implies 

l(l-w)L^$l^  ♦  0,  n  ♦  N. 

It  follows  that 

|<f,»>  -  li.  «Vv>| 

N-H» 

tends  to  zero,  uniformly  for  bounded  x  (cf.  definition  of  ♦). 

Proof  of  Theorem  2.  Let  f  e5  and  assume  that  supp  $  fl  15^  -  0.  If  the 
sequence  sn  e  converges  locally  uniformly  to  f  then  (2)  holds  and 
we  have 

<f . ♦>  -  llm  (ill,  l  <a(J)L(-l),  }>). 

»*■  K~ 
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Let  A  :■  +  (-gjr-)  denote  the  Laplace  operator.  Since  Ln 

together  with  all  its  derivatives  vanishes  at  infinity,  we  have 

<8n(j)Ln(*~J)’  ^  “  ^0^'  **n»  ^ 

m  j.  I  <  I  2  i”2— c/2  \  — J ij  •  / * \2+c/2/^*  i  \ x 


Xi 


Applying  Theorem  3  with  e  -  dist(supp  $,  1^)  >  0  we  have  for 
sufficiently  large  n', 

|<...>|  <  c(l  +  |j|2r2[sup(l  +  |j|2)C/2|8nO)|  Jci  +  ce)’n'l^l. 

It  follows  that  <£, ♦>  -  0. 


For  the  proof  of  Theorem  3  we  make  use  of  the  following  precise 

estimates  for  L  and  the  numbers  a  .  which  have  been  derived  In  [31. 
n  n,  j 


Theorem  4  [3,  Thm.  5.2], 


Lr(x)  -  ^(x)]  <  c(l  +  cdn(x))-n  . 


Proposition  1.  [3,  Prop.  5.2,  Lemmas  6.5,  6.6].  Set  J'  :» 


{*(!»!)»  *(2,-1),  i(-l,2)}.  Then,  for  (2ux)  e  «n,  we  have 


(l  +  c  dist(2nx,  r  fl  0n) J””' *  j  e  J 

(7)  la  (x) I  <  ■  (l  +  c  dist  (2*x,  —  j / 2 )  )~tl  ,  j  e  J' 

•  n,  J  J 

(l  +  c|j|  J-”',  j  e  Z2  \{J  U  j'  U  0}  . 

Proposition  2.  [3,  Prop.  5.1].  fln  depends  continuously  on  n  in  the 

Hausdorff  topology. 

The  reader  who  compares  these  statements  with  those  in  [3]  will  notice 
that  we  have  slightly  changed  the  notation.  Note  that  the  estimate  (6)  is 
stronger  than  the  assertion  of  Theorem  3  for  a  -  0  since  the  constants 
c  in  (6)  do  not  depend  on  the  distance  of  x  to 
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need  the  analogue  of  estimate  (7)  for  the  derivatives  of  a  . . 

n»  J 


Lean a  1.  For  any  5  >  0,  there  exist  constants  c^,  c  and  Oq(6)  such 


that  for  all  n'  >  nrt  and  2*x  e 
— — —  o  —  n 


[i  +  c  dist(2*x,  r  n  a  )]_n  ,  j  e 
n,j  n 


(8)  |Daan  .j(x)|  <  c6(l  +  cfi)n'  [l  +  c  diet(2*x,  -j/2)]_n‘,  j  e  J* 

.  U  +  c|j|  j  el2\{JUJ'  U0}. 


The  proof  of  this  Lemma  is  technical  and  we  postpone  it  until  the  end. 


Proposition  3.  Let  x'  ■  x  +  j  with  J  e  2  \  0  and  2*x  e  Then, 


for  any  6  >  0,  there  exist  constants  c^,  c  and  nQ(6)  such  that  for 


all  n’  > 


(9)  |Daan>J(x)|  <  cfi(l  +  c«)n'(l  +  cdn(2xx»))"°' 


For  a  ■  0  this  is  Proposition  5.4  in  [3).  There  we  bounded  the  terms  in 


square  brackets  on  the  right  hand  side  of  (8)  by  (1  +  cd  (2trx' ) )  which 


appears  on  the  right  hand  side  of  (9).  Clearly,  the  case  a  *  0  can  be 


treated  in  the  seam  way. 


Proof  of  Theorem  3.  Since 


Ln(2wx) 


Pn(2*x> 

-* -  *  l  2  a„  j(x>» 

Mn(2xx)  je*  n»J 


we  have,  for  a  -  1, 


(2x)Da£  (2ix)  -  -L  (2sx)2  l  D«a  (x), 
n  n  j#Q  nj 


For  arbitrary  o  *  0  if  follows  that 


(10) 


(2»)  a  D®L  (2wx) 


l 

J#0 


ei+S2<a 


*1. 

c-D  aL  (2*x)D 
p  n 


a-8.-B2 

\ (2wx)D  1  \  ,(x). 

a  n,j 


Let  us  first  assune  that  2wx  e  a  .  We  data  that  for  any  8  >  0 
there  exists  such  that 

(11)  |DaLn(2iix)|  <  cfi(l  +  c 8)n ' ( 1  +  cdn(2xx))"0,t  2»x  e  Qn  . 

i 

For  a  -  0  this  is  a  weaker  statement  than  the  assertion  of  Theorem  4 
Using  Induction  on  J  a|  it  is  sufficient  to  show  that 


J*0 


can  be  bounded  by  the  right  hand  side  of  (11).  Leona  1  yields 


|Daan ^(x)|  <  cfi(l  +  c«)n'  •  - 


(1  +  cdn(2»x))“  JtJUJ', 

(l  +  c|j|  J  e*2\  {J  U  J*  U  0) 


Sunning  this  Inequality  over  j  t  Ir  \  0  finishes  the  proof  of  (11). 


Secondly,  let  x*  ■  x  +  j  with  (2ix)  e  Then,  writing 

C  (2ix’)  -  C  (2wx)an  .(x), 
n  n  n,  j 

we  see  that 

Da£n(2ix»)  -  l  cftDBCn(2iix)D*"3an  .(x)  . 
n  *<<*  3  n  n»* 

Therefore,  by  (11)  and  Proposition  3,  D<1^n  can  be  estiaated  by 


1 


|DaLn(2wx,)|  <  c5(l  +c6)n'(l  +  cdn(2wx'))-n'  , 


x'  *  x  +  j ,  2wx  e  ft 


Theorem  3  easily  follows  from  the  estimates  (11)  and  (12):  Let  e  >  0 

and  assume  that  d^(2xx)  >  e.  We  choose  ng  so  that  dn(2vx)  >  d{](2xx)/2 

for  n'  >  ng.  Now  (11)  and  (12)  give  (4)  since  we  can  choose  6 

sufficiently  small. 


Proof  of  Lemma  1.  In  proving  (8)  we  make  use  of  ttye  symmetries  of  the 
mesh.  If  A  is  a  linear  transformation  which  leaves  the  set  J 
invariant,  we  have 


an,Aj(Ax)  "  an,j(x> 


where  n  is  the  appropriate  permutation  of  n.  Similarly, 


Aft  -  ft-, 
n  n 


From  this  one  can  check  (cf.  [3,  section  3])  that  one  may  assume,  by 
changing  n  if  necessary,  that  x  :»  (u,v)  lies  in  the  first  quadrant. 
Further,  since  the  roles  of  u  and  v  may  be  interchanged  and 
(2«x)  e  ftn  c  ft  ,  (cf .  Figure  3)  we  shall  assume  throughout  this  proof 


0  <  v  <  u  <  1/2  . 


By  definition  (1.5)  of  a  .  we  have 

n»  j 


■  V  .  *  *-  ^  •. 


|D°*n,j(u’v>l  i 


,"r9i 


,“2"62 


i"3'63 


n.+B,  «o+Bc  n 

|u+kl  1  4  jv+il  2  5  |u+v+k+il  3  6 


where  the  sum  is  taken  over  all  £  which  satisfy 


°  <  £v  <  ny,  v  -  1,2,3, 


1  *v  “  ai  +  «2 
\*1  1  z 


k  -  0 


"®3fv  in  ca8e  |  *  “  0 


k+i  *0  v  * 


This  last  restriction  comes  from  the  fact  that,  e.g.,  for  k  -  0,  the 


u  ^nl 


factor  1  is  equal  to  1,  hence  does  not  figure  in  the  differentia¬ 

tion.  To  estimate  the  individual  summands  in  (14)  we  consider  4  cases. 
Unless  (k,i)  c  {(0,-1),  (-1,1))  (cases  (ii)(b),  (c))  we  bound  each 
summand  (...)  on  the  right  hand  side  of  (14)  by 


cfi(l  +  c6)n'  max((  1  +  c|j|  p',  |«n )  • 


(i)  v  <  fi  <  1/4,  u  <  1/8: 


Using  the  inequality 


pjy  <  U  ♦  ceM  rl.  q  ez\0,  |p|  <  1/2  -  e  , 


we  obtain  the  estimate 


m 


"n2+e2 


-n3+e3 


|[— ]|  <c(l  +  c|k|)  l(l+c|i|)  (1  +  c|k+£|) 

,  r<0i+a4> 

with  c  involving  terms  like  u+k  ,  k  *  0,  which  are  bounded. 

Since  at  most  one  of  the  components  of  n  is  less  than  n',  this  implies 


(ii)  v  <  «  <  1/4,  u  >  1/8  : 


We  consider  several  subcases. 


(a)  (k,l)  *  {(0,1),  (-1,1)): 


We  have 


2+02  i  u  inli  u+v  in3 


This  can  be  estimated  as  before  unless  k  ■  -1  or  k+£  »  -1.  If  k  -  -1 
and  l  *  0,1  we  can  write  the  right  side  of  (18)  as 


.  l.M  2  P2|  ,  v  I  u+v-1  i"3 

c(*  +  c|  *|  )  l.+i-1+tl 


The  last  factor  is  less  than  (l  +  c|i|  )  and,  since  for  2*(u,v)  e  fln, 

|an  (-1  0)(tt’v)l  <  l*  ^  follows.  If  (k, i)  -  (-1,0)  it  is  easily 
seen  that  the  left  hand  side  of  (18)  can  be  bounded  by  c  |aQ  Q^(u,v) j . 
The  cases  k  +  1  ■  -1,  k  *  0,  -1  are  treated  similarly. 

(b)  (k.t)  -  (0,-1): 

Set  n^  :■  (nj,  n2“&2*  n3^*  Proposition  2,  there  exists 

nn  -  nn(a, 6)  so  that  the  boundaries  of  0  and  0  are  within  6  of 
u  u  n  Op 

the  boundary  of  the  limit  set  for  n'  J>  no  and  a1-1  sati8fyln8 


w 


(15).  Moreover,  y  n*  £nj}£2n,t  n'  >_  nQ. 


For  u+v  <  1/2  (and  n'  >  n^)  we  obtain,  using  also  Proposition  1, 


•lUc!*v<o.-i)<“*''>l 


<  c(l  +  c  dlst (2*x,  rn  (0t-i)  n  )) 

p  p 

<  c(l  +  c6)n'(l  +  c  dlst(2irx,  Tn  iQ  _x)  0  8j)‘ 


For  the  last  Inequality  we  have  used  that 
di,t(rn6,(0,-l)  0  flng’  rn,(0,-l)  n  i  6* 
For  1/2  <  u+v  <  1/2  +  fi  we  have 


v  in2“e2.  u+v  in3 


f  l  <  |  V  I  Z  Z|  u+v 

I1—1  -  cIt=v!  h=s=7 


v  In2"a2 


<4^p  -MiTSy"3. 


where  we  have  used  that  aQ  Q^(u,v)  <1.  By  our  assumptions  on  u 

and  v  the  alnlnun  can  be  estimated  by  (1  +  c6)n  .  Therefore,  if 

n.  >  cn'y  (8)  follows.  If  lim  n0/n '  -  0,  the  curve  P  ,n  , v  0  TT 
2  —  2  n,iu,-lj  n 

converges  to  the  segment  (2»(u,v):  u+v  ■  1/2,  (u,v)  >  (0,0)}.  Therefore 


we  nay  assume  that 


diet (2wx,  rn  (0  n  nj  <  c« 


for  n'  and  (8)  follows. 


(c)  (k,i)  -  (-1,1): 


We  have 


n2"*2 


-n, 


<  c(l  +  c)  ‘  (1  +  c( 1/2  -  u))  \ 


Since  T  .  ,  ,v  does  not  intersect  the  square  [0,*]  x  [-11,0],  this 
n, l-l, 

iaplies  (8). 

(ill)  6  <  v,  u  <  1/2  -  6:  Since  v  <  u,  we  have 


> 

(iv)  5  <  v,  u  >  1/2  -  5:  We  have 


u+v  i  p6 


u+v  |  p6 


Assume  e.g.,  that  n^  -  ain(n1,n2)*  Since  2 wx  e  Qn, 


K(-i.o»«l  -  <  1 


Froa  this  and  the  fact  that  u  >  1/2  -  6,  8  <  n  ,  we  have 

*  o  —  j 


H  a 

.  ,  n«  6  , 

|^r|6<|^|3  <  <■  +  cS>  • 


Coabining  the  above  estimates  yields 


|>  — 'I  ic«K.J<*>l(1  +cS)n' 
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